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ABSTRACT 
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I 

sites.  Cumulative  frequency  curves  suggested  that  there  were  two  major 
populations  for  all  metals  and  perhaps  a  third  and  smaller  one  for  Cd,  Cr, 
and  Mn.  Plots  of  metal  ^  Fe  indicated  no  anthropogenic  inputs  of  metals 
for  shelf  and  Chesapeake  Bay  channel  sites,  but  suggested  anthropogenic 
influences  for  all  metals  in  several  of  the  inshore  sites.  Enrichment 
factor  calculations  showed  enrichment  of  Cd,  Cu,  Pb,  and  Zn  with  respect 
to  average  crustal  abundances.  A  recommendation  of  this  study  for  evalu¬ 
ation  of  environmental  geochemical  metals  data  is  to  utilize  mean  concen¬ 
trations,  cumulative  frequency  plots,  and  metal  ^  Fe  and/or  enrichment 
factor  calculations  when  evaluating  the  pollution  status  of  sediments. 
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TRACE  ELEMENTS  AND  THEIR  ASSOCIATIONS 
IN  SEDIMENTS  OVER  THE  HAMPTON  ROADS  HARBOR 
AND  LOWER  CHESAPEAKE  BAY  AREA 


Introduction 

When  assessing  metal  sediment  concentrations  for  environmental  studies, 
one  major  problem  is  the  choice  of  methods  of  data  analyses.-  One  may 
attempt  to  evaluate  the  data  based  on  cibsolute  metal  concentrations,  or 
choose  between  a  variety  of  other  methods,  ranging  from  relatively  simple 
ones  like  elemental  ratios  to  more  sophisticated  ones  such  as  discrimi¬ 
nant  analysis. 

Helz  ^  al.  (1983)  analyzed  sediments  from  the  Chesapeake  Bay,  Baltimore 
Harbor,  and  Elizabeth  River  for  trace  elements  (Cr,  Mn,  Fe,  Co,  Ni,  Cu, 

Zn,  Cd,  and  Pb) .  They  calculated  enrichment  factors  (E.F.)  for  their  data 
using  Al  as  the  reference  metal  and  average  crustal  shale  values  from 
Turekian  and  Wedepohl  (1961)  by  the  relation: 

_  (X/Al)  sediment 

(X/Al)  Ecurth's  Crust 

where  X/Al  is  the  ratio  of  the  concentration  of  element  X  to  Al.  They  found, 
in  general,  enrichment  of  Cd,  Pb,  and  Zn  in  the  Bay  and  postulated  that  this 
resulted  from:  (1)  atmospheric  deposition,  or  (2)  unrepresentative  crustal 
abundance  data  -  perhaps  metal  concentrations  in  Chesapeaike  Bay  are  anomalous 
as  compared  to  the  average  crustal  abundances. 


Using  Fe  based  enrichment  factors,  Helz  e^  al.  (1983)  compared  their 
data  to  that  from  several  other  East  and  Gulf  Coast  estuaries  and  found 
enrichment  of  Cd,  Pb,  and  Zn  in  most  cases.  They  stated  that  the  Al 
concentrations  in  both  the  Hampton  Roads  and  Baltimore  Harbor  sediments 
were  high  when  ccxnpared  to  the  respective  main  Bay  areas.  In  Hampton 
Roads,  no  consistent  anthropogenic  effect  was  found  for  Cr  or  V  but  Zn 
was  highly  enriched. 

Goldberg  e^  a^.  (1978)  collected  and  analyzed  eight  cores  from  the 
Chesapeedce  Bay  (including  one  from  the  James  River) .  They  presented 
detailed  metal  concentrations  for  four  of  these  cores  and  concluded,  from 
the  eJasolute  concentrations,  that  surface  metal  contents  were  considerably 
enriched  over  values  below  50  cm.  This  trend  was  not  cleeurly  evident  in 
all  cases.  In  some  instances,  other  types  of  data  analysis  may  have 
provided  a  clearer  indication  of  whether  or  not  there  was  enrichment. 
Normalization  of  the  trace  metal  concentrations  to  either  Al  or  Fe 
may  have  been  very  useful  for  this  purpose. 

Sommer  and  Pyzilt  (1974)  used  metal  versus  aluminum/metal  ratios  as  an 
aid  to  evaluate  metal-sediment  associations  of  Chesapea)ce  Bay  sediments. 
They  determined  that  Al  was  associated  with  clay  minerals  in  Bay  sediment 
and  used  Al  as  an  indicator  of  clay  minerals.  Linear  relationships  were 
found  between  Cu  £uid  Al/Cu;  Pb  and  Al/Pb;  Cr  and  Al/Cr;  Mn  and  Al/Mn. 

Iron  did  not  show  a  linear  relationship.  Sommer  and  Pyzik  suggested  that 
because  of  the  linear  relationships  for  Al  and  Cu,  Pb,  Cr,  and  Mn, 
these  elements  were  associated  with  the  clay  mineral  portion  of  the 
sediment  and  that  an  alternate  distribution  mechanism  (sulfides) 


was  possible  for  Fe. 


Johnson  and  Villa,  Jr.  (1976)  determined  acid  extractable  metals  in 
Elizabeth  River  sediments.  They  found  that  concentrations  of  all  metals 
analyzed  (Cd,  Cu,  Cr,  Hg,  Pb,  Zn,  Al,  and  Fe)  were  two  to  ten  times  greater 
than  for  sediments  from  the  mid-Chesapeake  Bay.  Distribution  of  metals 
generally  reflected  heavy  industrial,  comnercial ,  and  domestic  inputs 
along  the  river.  Their  findings  are  reviewed  in  the  discussion  section. 

Johnson  and  Villa,  Jr.  (1976)  discussed  the  work  of  Sommer  amd  Pyzik 
(1974)  and  examined  their  data  for  metals  in  Elizabeth  River  sediments  by 
the  same  metal  v^  Al/metal  technique.  No  linear  relationships  were  found 
for  Fe,  Cr,  Pb,  or  Cu  with  Al.  They  then  concluded  that  either  Al  was  not 
associated  with  clay  minerals  in  the  Elizabeth  River  as  it  was  in  the  Bay 
sediments,  or  non-linear  relationships  are  indicative  of  anthropogenic 
sources  rather  than  natural  ones. 

It  seems  highly  unlikely  that  the  association  of  Al  with  clay  minerals 
(assuming  these  authors  mean  silicate  clay  minerals)  would  vary  as  postu¬ 
lated  cibove,  since  Al  is  an  inherent  and  abundant  component  of  the  clay 
minerals  structure.  Other  forms  of  Al  were  undoiibtedly  present  in  sediments 
since  the  extraction  techniques  used  by  these  workers  would  not  significantly 
dissolve  the  silicate  lattice.  Whether  or  not  these  other  forms  of  Al  would 
be  associated  with  silicate  clay  minerals  may  be  questionable,  but  it  is 
expected  that  these  forms  will  be  in  the  clay-sized  fraction  of  sediments. 
Most  of  the  Fe  in  sediments  will  be  in  the  forms  of  hydroxides  or  sulfides, 
depending  on  Eh,  and  in  the  clay-size  fraction.  The  most  relevant  con¬ 
clusion  as  to  the  association  and  distribution  of  Al  and  Fe,  with  respect 
to  clay  minerals,  is  that  the  forms  of  all  three  occur  predominantly 


in  the  clay-sized  sediment  fraction.  Many  sediment  forms  of  A1  and  Fe 
are  relatively  easily  extracted;  silicate  clays  are  not. 

Trefry  and  Presley  (1976a,  b)  used  metal  v£  Fe  plots  to  characterize 
unpolluted  sediments  and  to  identify  sediments  thought  to  be  polluted  from 
anthropogenic  sources.  They  did  not  use  metal  v£  metal/Fe  ratios  but 
simple  metal  v£  Fe  scatter  plots.  A  population  that  occurred  within  the 
^95%  prediction  interval  about  the  best  fit  linear  line  was  defined  as 
natural  or  unpolluted.  Metal  concentrations  occurring  outside  of  the  95% 
interval  were  postulated  to  indicate  anthropogenic  input.  Rule  (1982) 
discussed  sediment  metals  in  the  lower  Chesapeake  Bay  and  Hampton  Roads 
Harbor  area  as  evaluated  by  the  technique  of  Trefrey  and  Presley. 

Chester  and  Voutsinou  (1981)  used  dilute  HCl  to  extract  sediment  metals 
and  assessed  the  sediment's  pollution  status  by  examining  absolute  trace 
metal  concentrations.  They  suggest  that  spatial  contour  maps  of  surface 
metal  concentrations  are  a  simple  and  effective  manner  to  present  results. 
Although  proposing  the  use  of  this  simple  method,  they  point  out  that 
great  care  is  necessary  in  choosing  baseline  metal  levels. 

Absolute  amounts  of  trace  elements  are,  to  some  degree,  simply  a  function 
of  sediment  particle  size.  Therefore,  assessing  a  group  of  data  where 
sediment  particle  size  is  variable  roust  be  done  with  caution.  'Even  when 
the  grain  size  is  nearly  homogeneous  within  a  data  set,  the  choice  of 
baseline  data  must  be  done  with  care. 

Sinex  and  Helz  (1981)  noted  that  certain  un fractionated  seunples  from 
the  Chesapeake  Bay  contained  higher  trace  element  concentration  than  cor¬ 
responding  <63ym  fractions.  They  postulated  that  the  coarser  grained 
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fraction  contained  fecal  pellets,  heavy  minerals,  aggregates  of  fine¬ 
grained  materials,  or  Scind  grains  with  unusually  thick  coatings. 

Shrader,  Rule,  and  Furbish  (1977)  studied  the  distribution  of  several 
metals  in  different  size  fractions  of  stream  sediments .  They  found  that 
the  concentration  of  some  metals,  notably  Mn,  was  greater  in  the  sand 
and  clay  fractions  than  in  the  silt  sizes  with  the  highest  concentrations 
often  in  the  sand  sizes.  This  phenomenon  is  most  probably  attributable 
to  Mn-oxide  coatings  on  the  sand  grains. 

Generally,  the  sediment  fraction  that  contains  the  greatest  concentra¬ 
tions  of  trace  metal  is,  of  course,  the  finer  size.  A  very  desirable 
method  of  assessing  the  trace  metal  status  of  sediments  would  be  to 
analyze  the  <63ym  material.  In  this  manner,  the  variable  effect  of  grain 
size  will  be  eliminated.  Most  environmental  studies  utilize  an  unfrac¬ 
tionated  or  <2mm  portion  of  sediments  for  the  various  determinations. 

Analyses  of  the  <63ym  fraction,  while  desirable,  would  entail  additional 
time  and  labor  often  beyond  that  available  for  the  research.  The  most  useful 
technique  should  allow  use  of  chemical  data  from  unfractionated  samples,  be 
applicable  to  a  wide  range  of  sediment  types,  and  minimize  the  difficulty 
in  selection  of  baseline  or  reference  concentrations.  A  number  of  such 
techniques  do  exist;  this  paper  evaluates  some  of  them  utilizing  sediment 
metal  data  from  the  lower  Chesapeake  Bay,  Inner  Virginia  Shelf,  and 
Hampton  Roads  Harbor  System.  A  suggested  plan  for  the  data  analyses 
of  trace  metals  in  sediments  is  given. 
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Materials  and  Methods 

The  Port  of  Hampton  Roads,  Virginia,  is  located  within  a  major  metro¬ 
politan  area  encompassed  by  Norfolk,  Virginia  Beach,  Chesapeake,  Portsmouth, 
Newport  News,  and  Hampton  (Figure  1) .  Waters  of  this  harbor  are  along  one 
of  the  most  industrialized  coastal  areas  in  the  eastern  part  of  the 
United  States.  The  largest  military  port  in  the  world  is  located  here. 
Navigation  channels  of  this  port  system  are  maintained  by  the  U.S.  Army 
Corps  of  Engineers  (COE).  The  COE  dredges  an  annual  average  of  4.1  million 
m^  of  sediment  from  the  port  systems  channels .  Approximately  60%  of  this 
material  is  classified  as  mud,  clay,  and  silt  and  is  taken  primarily  from 
the  HamptonRoads/Elizabeth  River  portion.  These  are  the  most  active 
industrial  areas  in  the  port  system.  The  balance  of  the  dredged  materials 
consists  of  sand,  gravel,  and  shell  which  comes  from  the  Chesapeake  Bay 
Channels  (Pequegnat  et  al . ,  1978,  current  research). 

The  large  Naval  Base  is  located  adjacent  to  Sites  E  and  F;  Sites  G-I 
are  on  the  main  stem  of  the  Elizabeth  River  where  there  is  moderate 
industrial  activity  (fuel  facilities  and  chemical  companies) .  The  highly 
industrialized  area,  with  private  and  naval  shipyards,  oil  storage  and 
creosoting  facilities,  and  fertilizer  companies  occur  along  Sites  K,  L, 
and  M.  A  less  densely  developed  area  in  the  upper  reaches  of  the 
Elizabeth  River  are  represented  by  Sites  P,  Q,  S,  T. 

Sediments  from  the  highly  industrialized  segments  are  presently  disposed 
at  Craney  Island,  a  diked  disposal  area  (Figure  1) .  With  normal  maintenance, 
dredging  and  disposal  operations,  the  Craney  Island  site  will  be  available 
for  several  years  but  has  a  finite  lifetime.  Plans  are  also  being  made  to 


deepen  the  Port  of  Hampton  Roads  from  the  present  13.8m  depth  to  17m  to 
accommodate  deeper  draft  vessels.  Should  the  deepening  occur,  the  filling 
of  Craney  Island  would  be  greatly  accelerated.  A  multidisciplinary  study 
was  designed  to  investigate  if  open  ocean  disposal  of  sediments  from 
various  areas  within  the  port  system  would  be  ecologically  acceptable. 

Sediment  saunples  were  collected  to  use  in  various  types  of  bioassays, 
for  particle  size  analyses,  and  for  metals  analysis.  Sediment  samples 
collected  for  all  these  various  purposes  were  analyzed  for  trace  metals. 

One  hundred  and  seventy-nine  sediment  samples  were  collected  over  a 
three-year  period  from  forty-six  stations  in  the  navigation  channels  in 
the  lower  Chesapecike  Bay  (Thimble  Shoal  and  Cape  Henry  channels) ,  Hampton 
Roads  Harbor,  and  the  Elizabeth  River  (main  stem  and  Southern  Branch) . 
Samples  from  a  proposed  offshore  disposal  site  on  the  Inner  Virginian 
Shelf  about  27  km  from  the  Bay  mouth  were  also  taken  (included  in  the 
179  count) . 

The  seunpling  pattern  for  the  offshore  area  is  shown  in  Figure  1.  All 
other  samples,  in  the  navigation  channels,  were  teiken  at  one-mile  inter¬ 
vals  using  navigation  buoys  as  site  markers.  At  each  site,  four  samples 
were  taken  across  the  channel  between  the  buoys . 

Sediment  samples  were  collected  using  either  a  Shipek  grab  or  a  0.76m^ 
cleunshell  grab.  Care  was  taken  to  avoid  contamination  from  the  samplers 
and  the  sauries  were  placed  in  plastic  bags  for  transport  to  the  laboratory 
Sediments  were  dried  at  <40C°,  crushed  and  passed  through  a  2mm  stainless 
steel  sieve  and  stored  in  plastic  bags.  Sediment  grain-size  analysis  was 
performed  on  the  samples  collected  from  the  Chesapeake  Bay  channels  and  the 


offshore  site  as  a  part  of  the  contracted  study.  Analysis  of  grain-size 
of  the  Hampton  Roads  and  Elizabeth  River  scunples  were  conducted  where 
time  permitted.  Portions  (2-6g)  of  the  <2mm  seimples  were  extracted  as 
follows:  fifteen  ml  of  concentrated  (15.4  M)  HNO3  (redistilled)  were 

added  to  the  samples  and  then  heated  at  lOQOc  for  four  hours.  Five  ml  of 
30%  H2O2  were  added  and  the  mixtures  heated  for  an  additional  2  hours.  The 
extract  was  filtered  through  acid  washed  Whatman  No. 2  filter  paper  and 
diluted  with  ultrapure  deionized  water  to  a  final  acid  concentration  of 
10%.  The  metals  Cd,  Co,  Cr,  Cu,  Fe,  Mn,  Ni,  Pb,  and  Zn  were  determined 
by  flame  Atomic  Absorption  Spectrophotometry  using  a  Perkin-Elmer  603 
instrument.  All  data  were  calculated  on  a  dry  weight  basis.  Stastical 
calculations  were  conducted  using  the  Statistical  Packages  for  Social 
Sciences  (SPSS)  by  Nie  et  al.  (1975) . 


Results  and  Discussion 


Sediment  metal  concentrations  as  averaged  for  different  areas  are  shown 
in  Table  1.  The  areas  are  defined  by  site  groupings  in  Table  2.  Metal 
concentrations  are  generally  lowest  in  the  Norfolk  Disposal  Site  (NDS)  area. 
Except  for  Pb,  concentrations  of  metals  are  similar  for  the  NDS  area  and 
Outer  Thimble  Channel.  The  Inner  Thimble  Channel  sediments,  closer  to 
inshore  activities,  have  a  slight  but  consistent  increase  in  concentrations 
for  all  metals.  The  NDS  sediments  are  fine  sands;  the  Thimble  Channel 

sediments  range  from  coarse  sands  to  fine  material  which  is  responsible 

0 

for  seme  of  the  increase  in  sediment  metal  concentrations. 

Starting  in  the  Inner  Thimble  Shoal  Channel,  the  general  pattern  shows 
an  increase  in  metal  concentrations  toward  the  industrialized  portions  of 
the  Elizabeth  River.  Levels  of  Cd,  Cu,  Pb,  and  Zn  exceeded  average  crustal 
abundances  (Taylor,  1964,  Table  3)  in  several  sites  close  to  industrial 
activ'ities.  Concentrations  of  Cd  and  Pb  showed  the  greatest  elevations 
above  crustal  abundances  (>10X)  and  this  increase  occurred  for  a  greater 
number  of  sites  than  for  other  metals. 

Maximum  concentrations  of  Cd,  Cr,  Cu,  and  Pb  occurred  in  sediments 
from  the  most  highly  industrialized  area  in  the  mid-Elizabeth  River,  while 
maximum  values  for  Fe  and  Zn  were  found  in  the  lower  Elizabeth  River.  The 
highest  values  for  both  .Mn  and  Ni  occurred  in  the  inner  Hampton  Roads 
sediments.  The  concentration  pattern  (Table  1)  of  Mn  and  Ni,  and  to  a 
lesser  extent,  Mn  and  Co,  suggests  a  geochemical  association  between  these 
elem.ents  in  the  sediments.  Geochemical  association  of  Mn  with  Ni  and  Co  has 
frequently  been  reported  (Siegel,  1974;  Levinson,  1980). 
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Concentrations  of  all  metals  decrease  in  the  Uppermost  Elizabeth  River, 
reflecting  the  lower  level  of  industrialization,  but  values  in  these  sedi¬ 
ments  were  still  greater  than  those  in  the  Thimble  Shoal  Chcuinel  and  NDS 
area,  and  exceeded  the  average  crustal  abundance  for  Cd,  Cu,  Pb,  and  Zn. 

Although  there  were  no  statistically  significant  differences  in  average 
concentrations  of  trace  elements  extracted  from  Elizabeth  River  sediments 
during  the  present  study  (sampled  in  1979)  and  the  levels  measured  by 
Johnson  and  Villa,  Jr.,  1976  (scimpled  in  1974,  Table  4),  the  ranges  for 

I 

most  metals  for  the  present  study  were  much  more  narrow  than  for  that  of 
Johnson  and  Villa,  Jr.  A  more  severe  extraction  procedure  was  used  in  the 
current  study  and  the  higher  mean  values  may  have  resulted  from  the  removal 
of  a  greater  amount  of  metals  from  the  sediments.  Although  the  average 
Johnson  and  Villa,  Jr.'s  values  were  lower,  there  were  several  instances 
where  maximum-found  values  exceeded  those  for  the  present  study.  This  may 
have  resulted  from  slightly  different  sampling  patterns  or  may  indicate 
a  change  in  the  distribution  pattern  of  the  metals. 

There  is,  of  course,  a  variation  in  grain  size  (Table  2)  over  the  study 
area  which  causes  some  of  the  differences  in  trace  metal  levels;  however, 
variations  in  metal  concentrations  also  exist  in  similar  grain-sized 
samples  (such  as  in  the  Elizabeth  River) .  Such  observations  are  common  in 
geochemical  studies.  Analyzing  the  <63um  fraction  of  sediments  may  help  to 
resolve  this  problem.  Since  it  was  not  possible  in  the  current  study  to 
analyze  the  <63pm  fraction,  the  next  best  approach  was  to  use  a  method  of 
data  analysis  which  might  account  for  grain-size  effect.  If  a  data  assess¬ 
ment  technique  can  be  utilized  which  inherently  incorporates  the  influence 
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(effect)  of  grain-size  on  the  trace  metal  content,  then  there  would  be 
no  need  to  group  the  data  by  particle  size  and  a  larger  set  of  data 
can  be  compared. 

Cumulative  frequency  plots  of  metal  concentrations  have  long  been  used 
(Siegel,  1974)  to  determine  the  number  of  populations  present  in  a  group 
of  samples.  When  a  cumulative  frequency  curve  for  a  given  metal  (metal 
concentration  ^  number  of  occurrences  for  each  concentration)  is  con¬ 
structed  on  a  probcibility  plot,  a  single  population  is  defined  as  a  group 
of  points  which  form  a  straight  line.  Multiple  populations  within  a  single 
(given)  data  set  will  give  two  or  more  straight-line  segments  on  the  graph. 

Cumulative  frequency  plots  were  constructed  for  each  metal  analyzed 
utilizing  the  entire  data  set.  All  curves  had  two  distinct  straight-line 
portions  indicating  two  major  populations  present.  Plots  for  Cd,  Cu,  Cr, 
and  Mn  (Figures  2  and  3)  show  the  variety  of  curves  obtained.  The  popu¬ 
lation  defined  by  the  near-vertical  segment  (low  metal  concentrations) 
includes  the  Norfolk  Disposal  Site  and  Thimble  Shoal  Channel  sediments. 

The  nearly  horizontal  segment  (higher  range  of  concentrations)  includes 
the  Hampton  Roads  Harbor  and  Elizabeth  River  sediments.  The  portion  of 
the  curves  for  Cd,  Cr,  and  Mn  between  the  larger  near-vertical  and  near¬ 
horizontal  segments  may  indicate  that  a  third  population  is  present. 

The  data  were  then  divided  into  three  groups  based  on  geographic  areas, 
sedimentological  and  geochemical  characteristics,  and  cumulative  frequency 
curves.  The  three  groups  were  Norfolk  Disposal  Site  (OFFSHORE);  Thimble 
Shoal  Channel  (CHANNEL) ;  and  sites  from  the  Hampton  Roads  Harbor  and 
Elizabeth  River  (INSHORE) . 


L2 


It  is  assumed  that  for  natural  sediment  conditions  (no  anomalous  metal 
concentrations) ,  the  relationship  between  Fe  and  another  metal  will  form  a 
linear  trend  (Trefrey  and  Presley,  1976a,  b) .  Should  the  concentration  of 
Fe  change  because  of  changing  mineralogy,  sediment  particle  size  (which 
will  inherently  change  mineralogy)  or  other  naturally  occurring  phenomena, 
the  concentration  of  the  other  metal  will  change  with  a  constant  relation 
to  Fe.  If  a  scatterplot  of  metal  Fe  is  constructed  and  a  regression 
line  with  a  95%  prediction  band  is  plotted,  then  the  natural  geochemical 
population  of  a  given  data  set  is  defined.  That  means  there  is  a  95% 
probability  that  the  points  that  fall  within  the  prediction  interval 
belong  to  the  normal  population  and  those  that  are  outside  of  this 
band  are  from  a  different  (or  ancanalous)  population. 

Metal/Fe  plots  were  constructed  for  the  OFFSHORE,  CHANNEL,  and  INSHORE 
data  groups  and  also  for  all  data  for  each  metal.  Representative  plots 
are  shown  in  Figure  4. 

Even  though  the  metal  concentrations  vary  greatly  between  the  data  sets 
in  Figure  4,  the  same  scale  was  used  for  comparative  purposes.  Correct 
plotting  of  OFFSHORE  Cu  for  evaluation  as  a  single  data  set  is  shown  in 
Figure  5a.  The  slope  and  nature  of  the  curves  for  the  OFFSHORE  and  CHANNEL 
scunples  (Figure  4a, b)  are  very  similar  and  no  doubt  from  an  unpolluted 
environment.  Metal  concentrations  are  low  and  all  points  fall  within 
the  95%  prediction  interval. 

When  all  sites  of  two  data  sets  belong  to  the  Scune  population,  they  will 
fall  along  the  same  line  of  the  metal/Fe  curves.  Those  saunples  having 
coarser  sediments  with  low  metal  concentrations  will  lie  on  the  lower 
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portion  near  the  ordineuit;  those  samples  having  finer  sediments  with 
resulting  higher  metal  concentrations  will  lie  further  to  the  right  but 
along  the  same  line.  It  is  clear  (Table  1)  that  the  Thimble  Channel 
sediments  vary  in  metal  content.  Ihe  outer  channel  sites  are  generally 
more  similar  to  the  OFFSHORE  than  those  in  the  inner  channel.  However,  the 
OFFSHORE  and  CHANNEL  samples  form  one  population  for  Cu  and  if  plotted 
together  would  fall  along  the  same  line.  When  samples  are  from  different 
geochemical  populations  and  have  different  metal/Fe  ratios,  they  tend  to 
form  lines  of  different  slopes  when  plotted  on  the  same  graph.  The  Cr/Fe 
plot  in  Figure  5b  is  a  very  good  example  of  vrfiat  appears  to  be  two  popu¬ 
lations  in  one  data  set.  The  upper  group  of  samples  to  the  left  fits 
closely  with  the  points  for  a  similar  OFFSHORE  plot  (Figure  5c)  while  the 
other  points  lie  along  a  distinctly  different  curve.  The  OFFSHORE  and 
CHANNEL  data  for  Cr/Fe  horizontal  row  of  points  are  from  relatively  coarse, 
clean  sands  in  the  Thimble  Shoal  Channel;  they  have  very  low  and  similar 
Cr  concentrations. 

The  data  distribution  for  the  INSHORE  sediment  plots  (Figure  4c)  is 
totally  different  than  for  the  other  two  groups.  The  points  are  widely 
scattered  and  the  95%  prediction  interval  is  relatively  wide.  Even  so, 
several  points  are  above  this  band,  indicating  anomalously  high  values. 
Although  there  is  some  particle-size  variation  in  the  INSHORE  group, 
sediments  are  predominemtly  silt-sized  or  less.  The  anomalously  high 
values  are  from  sites  along  the  highly  industrialized  portion  of 


the  Elizeibeth  River, 


Each  of  the  metal/Fe  (Cu/Fe)  plots  in  Figure  4a,  b,  c  have  included 
similar  grain-sized  material.  The  coarser,  low  metal  content  OFFSHORE  and 
CHANNEL  sediments  gave  similar  curves  with  few  anomalous  points.  The 
INSHORE  data,  with  mostly  silt  and  clay-sized  materials  gave  a  scatterplot 
with  widely  spread  points,  a  wide  prediction  interval,  and  several 
anomalous  points. 

When  data  from  all  sites  are  plotted  on  the  saune  graph  (Figiure  4d)  ,  less 
satisfactory  results  are  obtained.  One  obvious  difference  is  that  the  95% 
prediction  band  has  been  compressed.  The  large  number  of  OFFSHORE  and 
CHANNEL  samples  that  plot  together  near  the  ordinant  along  the  x-axis  had 
the  effect  of  compressing  the  confidence  band.  Even  so,  only  three  addi¬ 
tional  scimples  fall  above  the  line  as  being  described  as  having  anomalously 
high  Cu  values. 

As  a  more  direct  test  of  the  influence  of  particle  size  on  the  metal/Fe 
plots,  seunples  for  \^ich  particle  size  analysis  had  been  conducted  were 
selected.  Figure  6a  shows  plots  of  Zn  Fe  for  all  data  available  and 
only  for  those  samples  having  >15%  clay  (Figure  6b) .  Because  of  the  nature 
of  the  study  area  sediments,  most  of  these  data  are  from  the  INSHORE  sites. 
The  INSHORE  data  are  shown  in  Figure  6c) .  The  plot  of  Zn  Fe  for  >15% 

clay  samples  gives  a  regression  line  and  95%  prediction  interval  which  fit 
the  data  much  better  than  for  the  plot  of  all  data.  Plots  for  other  metals 
when  grouped  according  to  clay  or  silt  and  clay  content  were  similar  to  that 
in  Figure  5.  The  more  narrow  confidence  band  is  again  evident  in  the  all 
data  plot.  It  appears  that  the  potential  difficulty  with  the  all  data 
plots  is  with  the  great  number  of  data  points  (from  OFFSHORE  and  CHANNEL) 


lying  together  near  the  ordinant  and  not  with  the  particle-size  difference. 
Improved  evaluation  of  the  lower  concentration  samples  (OFFSHORE  and  CHANNEL 
sands)  can  be  better  accomplished  by  plotting  them  separately  or  as  one 
group  (Figure  5b)  but  evaluation  of  the  INSHORE  data  set  was  affected 
only  slightly  by  the  all  data  plot. 

Another  method  for  assessing  geochemical  data  is  the  enrichment  factor, 

EF  (Zoller  et  al . ,  1974;  Sinex  and  Helz,  1981).  The  EF  is  defined  for 
this  data  analysis  as : 

I 

_  (X/Fe)  sediment _ 

(X/Fe)  Earth's  crust 

where  X/Fe  is  the  ratio  of  the  concentration  of  element  X  to  Fe.  Anthro¬ 
pogenic  inputs  of  Fe  are  relatively  small  as  compared  to  natural  sources 
(Helz,  1976)  so  Fe  was  chosen  as  the  element  for  normalization.  Aluminum 
is  also  used  as  the  element  for  normalization  (Zoller  et  al. ,  1974;  Helz 
et  al. ,  1983)  and  thought  by  these  workers  to  perhaps  be  more  desireable 
than  Fe.  Since  Al  data  were  not  available  for  all  sediments  in  the  current 
study,  Fe  and  Al  could  not  be  compared,  nor  Al  used.  For  the  samples 
vfhere  both  Fe  and  Al  analysis  were  available,  plots  of  Fe  ^  Al  gave 
a  correlation  of  0.90  (Figure  6d) . 

The  major  assumption,  or  perhaps  limitation,  of  this  method  is  the  use 
of  average  crustal  edsundance  data.  A  study  area  may  be  either  enriched 
or  depleted  in  relation  to  the  average  crustal  abundance  and  still  be 
free  of  anthropogenic  inputs.  This  method  does,  however,  give  a  common 
reference  point  for  comparison.  A  value  of  1  would  mean  neither  enrichment 


nor  depletion  in  relation  to  the  Earth's  crust.  The  crustal  abundance 


data  of  Taylor  (1964)  were  used  to  compute  the  EF  for  all  data  in  this 
paper.  Sinex  and  Helz  (1981)  state  that  the  EF  is  not  very  sensitive  to 
the  choice  of  reference  material  and  if  average  shale  data  were  used 
instead  of  average  crustal  abundance,  the  enrichment  factors  would  change 
by  less  than  10%  except  for  Cd,  Pb,  and  Zn  which  would  decrease  ~40-50%. 

Enrichment  factors  for  metals  from  the  study  area  sediments  are  shown 
in  Table  5.  The  enrichment  factors  were  computed  for  individual  samples 
within  each  location  group  and  then  individual  EFs  were  averaged  to 
obtain  the  group  value.  Elements  that  show  definite  enrichment  are  Cd, 

Cu,  Pb,  and  Zn.  The  greatest  enrichment  occurs,  as  expected,  in  the 
heavily  industrialized  Elizabeth  River.  Cadmium  and  Pb  show  the  greatest 
enrichment  with  Zn  next  and  the  least  for  Cu.  The  enrichment  patterns  for 
these  four  elements  are  similar;  almost  identical  for  Cu,  Pb,  and  Zn,  and 
slightly  different  for  Cd.  The  seeningly  high  enrichment  factors  for  Cd 
and  Pb  for  the  OFFSHORE  and  Thimble  channel  sediments  are  probably  because 
of  the  relatively  low  concentrations  for  Fe.  Cadmium  and  Pb  concentrations 
in  these  sediments  were  very  low;  several  will  have  to  be  reanalyzed  with 
larger  aliquots  to  give  values  above  detection  limits.  Given  the  nature 
of  the  sddiments  (sands)  emd  the  very  low  Fe  values,  these  enrichment 
factors  should  be  cautiously  interpreted.  Even  though  the  OFFSHORE  s^unples 
are  defined  as  being  enriched,  the  absolute  concentrations  of  Cd  and  Pb 
obviously  indicate  no  need  for  environmental  concern.  The  enrichment 
factors  for  Cd,  Cu,  Pb,  and  Zn  are  within  the  range  of  values  computed  by 
Sinex  and  Helz  (1981)  for  various  East  and  Gulf  Coast  estuaries.  Enrichment 
factors  for  Co,  Cr,  Mn,  and  Ni  were  at  or  below  one. 
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Mean  enrichment  factors  from  the  data  comparison  with  Johnson  and  Villa,  Jr 
(1976)  are  presented  in  Teible  6.  Strong  enrichment  of  Cd,  Pb,  and  Zn  is 
evident  in  both  sets  of  data.  However,  there  was  an  apparent  decrease  in 
the  enrichment  factor  for  Cd  from  1974  to  1979  in  the  lower  portion  (Main 
Stem) .  This  decrease  most  li)cely  resulted  from  reduced  input  of  wastes 
and  wastewater  from  a  chemical  plant  along  the  river.  Factors  for  Cr  are 
unity  (or  near)  in  all  cases  and  there  has  been  no  change  for  Cu,  Pb, 
or  Zn  since  1974. 

Summary  and  Conclusions 

The  mean  concentrations  for  extractable  sediment  metals  were  evaluated 
with  respect  to  average  crustal  abundances.  Levels  of  Cd,  Cu,  Pb,  and  Zn 
exceeded  average  crustal  abundances  for  most  of  the  study  sites.  Mean 
values  for  Cd,  Cr,  Cu,  Fe,  Pb,  and  Zn  were  not  statistically  different 
from  scunples  measured  in  1974  in  the  Elizabeth  River.  The  range  of 
values  for  the  current  study  was  much  more  narrow  than  in  1974. 

Cumulative  frequency  curves  suggested  that  there  were  two  major  popu¬ 
lations  for  all  metals  and  perhaps  a  third  and  smaller  one  for  Cd,  Cr, 
and  Mn. 

Plots  of  metal  v£  Fe  indicated  that  sediments  from  the  NDS  site  and 
The  Thimble  Shoal  Channel  were  not  anthropogenically  affected  and,  for 
most  metals,  were  slightly  different  in  their  geochemical  nature.  These 
plots  for  INSHORE  data  were  very  different  and  indicated  that  several  sites 


were  enriched  well  cdaove  natural  levels.  Grouping  of  sediments  as  a  function 


of  particle  size  gave  less  scatter  of  the  points  and  identified  anomalously 
high  samples.  Based  on  the  present  results,  segregation  of  samples  based 
on  general  particle  size  (i_.e^-  #  sand  or  silt  and  clay)  before  constructing 
metal  Fe  plots  is  recommended. 

When  the  data  were  evaluated  using  the  enrichment  factor  (EF) ,  Cd,  Cu 
Pb,  and  Zn  showed  enrichment  with  respect  to  crustal  abundances.  Because 
of  the  very  low  Cd  and  Fe  values  for  the  NDS  area  which  resulted  in  high 
EFs,  it  is  recommended  that  both  sample  concentrations  and  EF  data  be 
presented  in  geochemical  studies. 

Both  the  metal  Fe  plots  and  the  EF  identified  anomalous  metal  concen¬ 
trations.  The  metal  Fe  plots  identified  anomalous  concentrations  for 
a  greater  number  of  elements  than  did  the  EF.  One  advantage  of  the  metal 
vs  Fe  plots  is  that  they  do  not  rely  on  average  2d3undance  data  which  may 
not  be  appropriate  for  the  study  area  sediment. 

Metal  concentrations,  metal  V£  Fe  plots  and  enrichment  factors  all 
identified  +-he  saune  sections  of  the  study  area  as  most  affected  by 
anthropogenic  sources.  Each  of  these  means  of  data  evaluation  contributes 
unique  and  useful  information  and  all  should  be  included  in  data  evaluation 

A  suggested  approach  to  geochemical  data  analysis  is  to:  (1)  evaluate 
metal  concentrations  using  appropriate  crustal  abundances;  (2)  construct 
cumulative  frequency  plots;  (3)  construct  metal  Fe  plots  utilizing 
particle  size  grouping;  and/or  (4)  calculate  enrichment  factors  choosing 
appropriate  reference  data.  In  all  instances,  the  metal  concentrations 
as  well  as  plots  or  enrichment  factors  should  be  presented. 


Metal  concentration  values  from  this  and  other  studies  (such  as  Sinex 
and  Helz,  1981)  are  obtained  from  partial  sediment  extractions  while  the 
crustal  cdsundances  reflect  total  values  for  a  given  rock  material.  It 
should  be  noted  that  the  EFs  so  calculated  are  conservative  estimates 
and  may  actually  be  greater  than  the  calculated  values.  Interpretation 
of  the  EFs  must,  therefore,  be  made  with  the  type  of  extraction  in  mind 
and  may  vary  from  study  to  study.  Interpretation  of  the  results  based 
on  the  extraction  method  can  lead  to  estimation  of  the  potential 
environmental  impact  based  on  enrichment  of  a  given  metal. 
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TABLE  2.  Site  legend  and  description  of  sediment  type. 


AREA 

SITES 

INCLUDED 

#  OF  SEDIMENT 

SAMPLES  TYPE 

NORFOLK  DISPOSAL  SITE 

Numbered  Sites 

49 

fine  sands 

THIMBLE  CHANNEL-OUTER 

AA  thru  EE,U 

35 

coarse  sands 

THIMBLE  CHANNEL- INNER 

FF  thru  II, V, A 

26 

coarse/fine  sands; silts 

NEWPORT  NEWS  CHANNEL 

B,C,D 

13 

silts/clays; sands ; shell 

INNER  HAMPTON  ROADS 

E,F 

9 

silts/clays 

LOWER  ELIZABETH  RIVER  (MAIN) 

G,H,I 

12 

clays/silts 

MID-ELIZABETH  RIVER 

K,L,M 

12 

clays/silts 

UPPER  ELIZABETH  RIVER 


P,C.S,T 


16 


clays/silts 


•rV 


TABLE  4.  Means  and  standard  deviations  of  trace  metals  in 
Elizabeth  River  sediments  from  the  current  study 
and  from  Johnson  and  Villa,  Jr.  (1976)*. 


MAIN  STEM 


METAL  (mg/kg)  1974,  EPA^  1979,  RULE^ 


SOUTHERN  BRANCH 


1974,  EPA^  1979,  RULET 


5. 3+5. 5 


47+26 


38+43 


70+68 


386+404 


2.0+J.l 

2 . 3+1 . 6 

2. 3+1.6 

52+6.6 

40+30 

57+17 

75+36 

82+91 

130+69 

36208+5512 

25092+8145 

31825+^9736 

93+39 

106+101 

149+62 

459+197 

287+263 

363+197 

Sites  were  chosen  from  the  current  study  to  correspond  with  the  areas 
sampled  by  Johnson  and  Villa,  Jr.  (1976). 

Sampling  for  the  Johnson  and  Villa,  Jr.  study  was  conducted  in  1974, 
and  in  the  Elized^eth  River  for  the  current  study  in  1979. 


/  >> 


n  IQ  *0 

0)  c 

o  u  *9 

o  Vt 

iP  c  ^  , 

0  +  I 


TABLE  6.  Comparison  of  enrichment  factors  from  the  current 
data  with  that  from  Johnson  &  Villa,  Jr.  (1976) 
for  the  Elizabeth  River. 


MAIN 

STEM 

SOUTHERN 

BRANCH 

METAL 

1974,  EPA 

1979,  RULE 

1974,  EPA 

1979,  RULE 

59  +741 

15  +7.9 

23  +12 

19  +9.4 

1.0+  0.9 

0.8+0. 1 

0.8+  0.4  ■ 

1 . 0+0 . 1 

1.4+  1.5 

2. 3+1.4 

3.1+  2.7 

4 . 0+1 . 3 

12  +12 

12  +4.2 

18  +14 

21  +4.0 

11  +11 

10  +3.2 

9+6 

8  +2.9 

Ranges  a^^e  Standard  Deviation. 


FIGURE  1.  Study  Area 


INSHORE  SEDIMENT  EXTR9CTI0NS  C0NF=95/: 
CORRELRTION  COEFFICIENT  =  .6198 
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C3  total  ORTfi  SEDIMENT  EXTRRCT'ONS  CC'-- 

r*  CORRELATION  COEFFICIENT  =  .7584 

*  3 


FIGURE  4.  concluded 


INSHORE  SEOIHENT  EXTRACTIONS  COtirrtS-. 
CORRELATION  COEFFICIENT  =  .8035 


*1 


total  OATR  SEOIHENT  EXTRACTIONS  CONFrgsz 
CORRELATION  COEFFICIENT  =  .9030 


FIGURE  6.  concluded 


